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CN ' Abstract 

f~| , We consider a Ginzburg-Landau partial differential equation in a bounded inter- 

' val, perturbed by weak spatio-temporal noise. As the interval length increases, a 

' ' , transition between activation regimes occurs, in which the classical Kramers rate di- 

verges |MS01| . We determine a corrected Kramers formula at the transition point, 
yielding a finite, though noise-dependent rate prefactor, confirming a conjecture by 
' Maier and Stein [MS03| . For both periodic and Neumann boundary conditions, we 

p I , obtain explicit expressions for the prefactor in terms of Bessel and error functions. 

c3 . 1 Introduction 

S: 

'— Weak noise acting on spatially extended systems can cause a wide range of interesting 

I phenomena. In particular, it can induce rare transitions between states which would be 

J> ' otherwise invariant, e.g. nucleation of one phase within another |Lan67j . micromagnetic 
domain reversal jNee491 IBra931 IBNROO| , pattern nucleation in electroconvection |CH93] , 

\^ , instabilities in metallic nanowires |BSS05j . and many others. The rate of such transitions 

CN I for weak noise intensity e is in general governed by Kramers' lav^ T ~ Tq exp{— AM/^/e}, 

Q^ ■ where the activation energy AW is the energy difference between stable and transition 

^ , states, and the rate prefactor Fq is related to second derivatives of the system's energy 

^ I functional at these states |Eyr35 , IKra40j . 



(N 



In a series of recent works |MS01t IMSOSj ISte04j , Maier and Stein studied transition 
rates in a Ginzburg-Landau partial differential equation on a finite interval, perturbed 
I by space-time white noise. They discovered the striking fact that as the interval length 

' approaches a critical value, which depends on the boundary conditions (b.c), the rate 

prefactor Tq diverges. Although this divergence is reminiscent of the behavior of certain 
thermodynamic quantities at phase transitions, it has a different origin |Ste05j : It is due 
to the fact that the Kramers law only takes into account the effect of quadratic terms in 
the energy functional on thermal fluctuations, while at the critical length some quadratic 
terms vanish due to a bifurcation, and higher-order terms come into play. 

Maier and Stein conjectured [MSOlj that the actual rate prefactor at the bifurcation 
point behaves like Fq — Ce~", for some constants C,a > 0. Until recently, Kramers rate 
theory was not sufficiently sharp to allow for the computation of these constants. Based 
on a new approach by Bovier et al [BEGK04] . we developed a method allowing to compute 
rate prefactors for potentials with nonquadratic transition states [BGOSbj . The aim of this 
Letter is to illustrate the method by determining the constants C and a in the case of the 
Ginzburg-Landau equation. 

^Throughout this Letter, the notation a b indicates that lim^^o o-/b = 1. 
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2 Model 



Consider a one-dimensional classical field (l){x,t), subjected to the quartic double-well 
potential energy function 

V{<P) = \<P' - \e , (2.1) 

to diffusion and to weak space-time white noise. Its evolution is given by the stochastic 
partial differential equation (SPDE) 

dt(t>{x, t) = d^^(bix, t) + cj){x, t) - cj){x, tf + V2^ C{x, t) , (2.2) 

where ^{x,t) denotes space-time Gaussian white noise, i.e., formally, 

H^{xiM)i{x2M)} = - X2)6{ti - ta) . (2.3) 

Here we consider the case of a bounded interval x £ [0,L], and either periodic or Neumann 
b.c. with zero flux, i.e., dx(p{0,t) = dx4>{L^t) = 0. 

Note that ^(x, t) can be rigorously deflned by independent white noises acting on each 
Fourier mode. For periodic b.c, this leads to setting 

^ oo 

cf>{x,t) = -^ E '^fc(i)e'"''"/^ (2.4) 

k=—oo 

and substituting in (12. 2p . The resulting system of stochastic differential equations (SDEs) 
is given by 

^k = ->^k4>k- 2 E (t)ki4'k24>k3 + V2eW[''\ (2.5) 

^1,^2,^3: ki+k2+k3=k 

where = — 1 -|- (27rfe/L)^, and the W^''^ are by deflnition independent Wiener processes 
(see for instance |Jet 86] for a discussion of the equivalence of different approaches to 
SPDEs). In the case of Neumann b.c, setting 

1 , . . A2" °° 



(^{x, t) = -j=Mt) + \ tYI ^k{t) cos{7Tkx/L) . (2.6) 



k=l 

yields a similar system of SDEs. 

With the SPDE (j2.2p we associate the energy functional 







i(0'(x))2 + y(0(x)) 



dx . (2.7) 



For both periodic and Neumann b.c, the uniform configurations (j)± = ±1 are stable 
stationary configurations of the system without noise. Both are minima of the energy 
functional, of energy 'H[(p±] = —L/A. In terms of the Fourier coefficients, for periodic b.c, 
the potential energy is given by 

1 °° 1 

WM=W[{(/>fc}] = - ^ Afc|0fc|2 + — Mk^Mk,, (2.8) 

k=—OD ki+k2+k:j+ki=0 

and a similar relation can be obtained for Neumann b.c. 

The value of the activation energy AW for this model is well known [FJL821 IMSOl] . 
The Kramers rate prefactor Fq, however, has only been determined for parameters L in 
certain ranges, excluding bifurcation values of the model [MS03| . 
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3 Transition states and activation energy 



The activation energy is the potential energy difference between the initial stable state 
and the transition state (pt . The latter is defined as the configuration of highest energy one 
cannot avoid reaching, when continuously deforming (p- to while keeping the energy as 
low as possible. The transition state is a stationary state of the energy functional, that is, 
it satisfies (j)"{x) = — (/>t(x) + (j)t{x)^. In addition, the Hessian operator S'^TL/Scj)^ must have 
a single negative eigenvalue at (j)^. The corresponding eigenfunction specifies the direction 
in which the most probable transition path approaches the transition state. 

The shape of (j)t depends on whether the bifurcation parameter L is smaller or larger 
than a critical value, the latter depending on the chosen b.c. [MSOlj . 



Periodic b.c. For L ^ 27r, the transition state is the identically zero function, which 
has energy zero. The activation barrier has thus value AVF = L/4. 

For L > 2tt, there is a continuous one-parameter family of transition states, of so-called 
instanton shape, given in terms of Jacobi's elliptic sine by 



2m 



■ sn 



+ ip,m 



m + l \y/m + 1 

Here if is an arbitrary phase shift, and m G [0, 1] is a parameter related to L by 

4Vm + 1 K(m) = L , 



(3.1) 



(3.2) 



where K(m) denotes the complete elliptic integral of the first kind. Note that m — > 0"^ as 
L approaches the critical length 2tt from above. Computing the energy of any instanton 
transition state (j3.ip . one gets |MS01| the activation barrier 



AW = n\ 



anstj 



1 



3^/r 



m 



8E{m) 



[1 — m){3m + 5) 
1 + m 



K(m) 



(3.3) 



where E(m) denotes the complete elliptic integral of the second kind. 



Neumann b.c. In this case, the identically zero solution forms the transition state for 
all L ^ vr, so that the activation barrier has again value AW = L/4. 

For L > vr, there are two transition states of instanton shape, given by 

/ 2771 I X \ 

4>nvst,±{x) = ±\ TT™ / , -I +K(m),m , (3.4) 

V m + 1 V ^771 +1 / 

where the parameter m G [0, 1] is now related to L by 

2Vm + 1 K(m) = L . (3.5) 

In this case we have m ^ 0"^ as L approaches the critical length vr from above. The 
activation energy is simply half the activation energy (j3.3p of the periodic case [MSOlj . 



3 



4 Rate prefactor 



The rate prefactor Tq is usually computed by Kramers' formula |Eyr35 , rKra40| 

1 



det Ae 



det At 



|At,o| • (4.1) 



Here As = d'^TL/dcf^lcp-] denotes the linearized evolution operator at the stable state (/>_, 
At = d'^Ti./d(j)'^[4>t] denotes the linearized evolution operator at the transition state cpt, and 
At,o denotes the single negative eigenvalue of At. 

For instance, for Neumann b.c. and L ^ vr, the eigenvalues of At = — d^/ dx^ + 1 are 
given by = — 1 + (vr/c/L)^, k = 0,1,2 ... , while the eigenvalues of Ag = — d^/ dx^ — 2 
are of the form ry^ = 2 + {■Kk/L)'^. It follows |MS03| that the rate prefactor is given by 



1 

2^' 



^ 2+{7Tk/Ly _ 1 / sinh(^/2L) 

\ J--'- |-1 + (^A;/L)2| 23/47rV sinL ' ^'^ 
\ fc=0 



The striking point is that this prefactor diverges, like (vr — L)^^/^, as L approaches the 
critical value tt. In fact, this is due to the Kramers formula ()4.1I) not being valid in cases of 
vanishing det At. To confirm Maier and Stein's conjecture that the rate prefactor at L = vr 
behaves like Ce~" and determine the constants C and a, we have to derive a corrected 
Kramers formula valid in such cases. This can be done |BG08b| by extending a technique 
initially developed by Bovier et al |BEGK04] . which we outline now. 

Potential theory. For simplicity, consider first the case of d-dimensional Brownian 
motion , starting in a point x G M°'. Given a set ^ C M'^, the expected value wa{x) = 
E[Tj|] of the first time the Brownian path hits A is known |Dyn65| to satisfy the 
boundary value problem 

/^wa{x) = 1 X e A" , 

wa{x) = xeA. (4.3) 

The solution can be written as 

wa{x) = / GA-{x,y) dy , (4.4) 

where Ga'^ denotes the associated Green's function, satisfying IS.,j;GA'={x,y) = 5{x — y) and 
the b.c. (for instance, Gj^3{x,y) = 1/{Att\\x — y\\)). Similarly, let hA,Bix) = IP{t^ < t^} 
denote the probability that the Brownian path starting in x hits the set A before hitting 
the set B. It satisfies the boundary value problem 

AhA,B{x)=0 x£{AuBy, 
hA,B{x) = 1 X e A , 

hA,B{x) =0 x&B . (4.5) 

This, however, is also the equation satisfied by the electric potential of a capacitor, with 
conductors A and B at respective potential 1 and 0. If pa,b{x) denotes the surface charge 
density on the two conductors, we can write 

hA,B{x) = I GB<^{x,y)pA,B{y) <^y ■ (4.6) 

JdA 
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The capacity of the capacitor is simply the total charge accumulated on one conductor, 
divided by the potential difference, which equals one: 



cap^(S) = / pA,B{y) dy ■ (4.7) 
JdA 

The key observation is the following. Let C = Bs{x) be a ball of radius e around x, 
and consider the integral jQ^WA{z)pc^Aiz) dz. On one hand, using the expression (j4.4p 
of WA, symmetry of the Green's function and then (14. 6[) . one sees that this integral is 
equal to J^^ hc,A{v) dy. On the other hand, as wa does not vary much on the small 
ball C |BEGK04] . we can replace wa{z) by wa{x)., and the remaining integral is just the 
capacity. This yields the relation 

/ his^{x),A{^) 

= E[t%] = wa{x) ^ ^-^ -J— . (4.8) 

The interest of this relation lies in the fact that capacities can be estimated by a variational 
principle. Indeed, the capacity for unit potential difference is equal to the total energy of 
the electric field, 

cap^(5) = / ||V/iA,B(x)f dx = inf / ||V/i(x)f dx , (4.9) 

J(A\JBY J(A\JBY 

where the infimum is taken over all twice differentiable functions satisfying the b.c. in (14. 5p . 

If, instead of Brownian motion, we consider the solution of a d-dimensional SDE 
X = —'VTl{x) + \/2eWt, the above steps can be repeated, provided we replace A by 
the generator eA — VH • V of the equation (the generator is the adjoint of the operator 
appearing in the Fokker-Planck equation). The above relations remain valid, only with 
the Lebesgue measure replaced by the invariant measure e~'^(^y^ dx. Thus we have 

r = E[r|]-i ~ — '^'^'^B.ixM) ^ .^^^ 

%W,A(^)e-^(^)/^ dz 



IA<= 

where the capacity can be computed via the Dirichlet form 

cavA(B) = inf e / ||V/i(x)f e'^^^)/" dx . (4.11) 
J{AyjBY 

The denominator in (j4.10p can be easily estimated by saddle-point methods, using the 
fact that /iB^(x),A is essentially 1 in the basin of attraction of x and in the basin of 
A. It is equal to leading order to (27re)'^/^ q-'H{x)/£ /^det((52'H/(^x2)(x). A good upper 
bound of the denominator in ()4.10p is obtained by inserting a sufficiently good guess for 
the potential h in (j4.1ip . Assume, e.g., that near a transition state at 0, the energy has 
the expansion 

1 1 "^'^ 

n{x) = --\Xo\xl + n(xi) + 2 5Z ^i^i + ' • • ' ^^'l^) 

where u(xi) corresponds to the possibly neutral direction in which a bifurcation occurs. 
Choosing /i(x) = /(xq) where e/"(xo) — dxQ'H{xo,0, . . . 0)/'(xo) with appropriate b.c. and 
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Figure 1. Rate prefactor Tq as a function of L/n for Neumann b.c. and different values 
of noise intensity s. 



substituting in (j4.1ip yields 



1 /(27r£)'^-i|Ao 



2tt 



A2 . . . Xd-i J- 



-n{xi)/e ^ 



(4.13) 



A matching lower bound for the capacity can be obtained by a slightly more elaborate 
argument, see [BGOSbj for details. If u{xi) = ^Xixf, the integral has value y^27re/Ai and 
we recover the usual Kramers formula. However, (j4.13p applies to other cases as well, e.g. 
a quartic u{xi). 

We now return to the SPDE (j2.2p . We apply the above theory first to a finite- 
dimensional approximation of the system (j2.2p . obtained either by truncation of high 
wave numbers in its Fourier transform (|2.5p . or by replacing the system by a discrete 
chain [BFGOTal IBFG07b , and then taking the limit. A difficulty is that the error terms 
will depend on the number of retained modes (see, for instance, |Liu03j for estimates on 
the convergence rate of the spectral approximations). Thus the results below are for now 
only formal. The error terms in the capacity can, however, be controlled |BG08a| . 



Neumann b.c. The potential energy along the normalized eigenvector in the bifurcating 
direction vi{x) = ^/2cos{t:x / L) is 



u{4)i) = H[(t)ivi] = L 



(4.14) 



Evaluating the integral in (|4.13p . we find |BG08bj that for L ^ tt the corrected Kramers 
prefactor to leading order is given by 



1 



Ai 



Ai 



h{V2L) 



' sm. 



sinL 



(4.15) 



where Ai = —1 + (vr/L)^ and is a universal scaling function, given in terms of the 
modified Bessel function of the second kind by 



a(l + a) 
8^^ 



(4.16) 
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For L <^ IT, since \I'_|_(a) tends to 1 as a —S' oo, we recover to leading order the rate (|4.2p . 
For TT — L of order -y/e, however, the correction terms come into play, and the factor \/Ai 
in the numerator of (|4.15p counteracts the divergence of the prefactor ()4.2p . In particular, 
we have 

r(i/4) 



^hm Fo c -^A_^^sinh(^vr) e-'/' . (4.17) 

For L > IT, the rate prefactor is harder to compute, because the transition states are not 
uniform. The computation can nevertheless be done |MS03| with the help of a method 
due to Gel'fand, with the result 



1 



sm. 



vr ' 



h{V2L) 



V2\{l-m) K(m) - (1 + m) E(m) 



(4.18) 



1 



where no 
L by ([33]). As L 

(L — 7r)~^/^. Proceeding as above, we find [BG08b] that the corrected prefactor is obtained 



— m + 1 is the negative eigenvalue of At , and m is related to 
7r+ (that is, m — > 0^), this expression again diverges, namely like 



by multiplying (|4.18p by 



2 V /"I + \/3e/4L V\/3e/4L 



(4.19) 



Here is again a universal scaling function, given in terms of modified Bessel functions 
of the first kind I±i/4 by 



*_(a) 



7ra(l+_a) 2/64 
32 



a- 



I \ 



a2\ 



(4.20) 



which converges to 2 as a — > oo, and /Ui is the second eigenvalue of At. We can in fact 
avoid the computation of this eigenvalue. Indeed, near the bifurcation a local analysis 
shows that ni = — 2Ai + 0(Af) = 3m + 0{m?)^ while further away from the bifurcation, 
the quotient in (j4.19p is close to 1. One can thus replace //i by 3m in (j4.19p . only causing 
a multiplicative error 1 + ©(e^/^). The resulting behavior of the prefactor Fq as L crosses 
the critical value vr is shown in Fig. [TJ 



Periodic b.c. For L ^ 27r, the transition state is uniform, and the computations are 
analogous to those in the previous case. The eigenvalues at the stable and transition states 
are now given by = — 1 + (27rfe/L)^ and % = 2 + (27rA;/L)^ with k e Z , and are thus 
double except for A; = 0. This implies that the integral in (14.13P is to be replaced by a 
double integral over the subspace of the two bifurcating modes [BGOSaJ . The result is 



1 



Ai 



Ai 



sinh(L/\/2) 
2tt X-^ + ,/3e/AL^^\^/3e/ALj sin(L/2) 



(4.21) 



where the scaling function is now given in terms of the error function by 



^+(a) 



''(l + a)e"'/8 



[1 + erf(-2-3/2^)] 



(4.22) 
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As ^I'^ converges to 1 as a ^ oo, for 27r — L » -^e, we recover the usual Kramers prefactor, 
which diverges as (27r — L)~^ as L ^ 2tt^ . However, as L approaches 27r, the correction 
terms come into play and we get 

^ sinh(-v/27r) 1/9 . 
lim To ~ )^ e~^'^ . 4.23 

For L > 27r, we again have to deal with a non-uniform transition state (/>t. An additional 
difficulty stems from the fact that transition states form a continuous family, so that the 
Hessian at always admits one vanishing eigenvalue. This eigenvalue can be removed 
by a regularization procedure due to McKane and Tarlie |MT95] . which has been applied 



in the case of an asymmetric potential in |SteQ4j . The computations are similar in the 
symmetric case |Ste] . and yield a rate prefactor per unit length 



To l^ol / 2m(l — m) sinh^(L/-v/2) _yi2 



L (27r)3/2y (l+m)5/2|K(m)-i±^E(m) 



e-'^' , (4.24) 



with Ayjm + 1 K(m) = L and the same /iq as for Neumann b.c. The factor e"-'^/^ reflects the 
fact that nucleation can occur anywhere in space |Ste04| . The prefactor now converges to 
a finite limit as L ^ 27r"^, which differs, however, by a factor 2 from ()4.23p . This apparent 
discrepancy is solved by applying the corrected Kramers formula, which shows that ()4.24p 
has to be multiplied by a factor 

/ 3m \ 

^ hr^^ (4.25) 



where ^>(x) = ^[1 + erf(x/\/2)]- The resulting rate prefactor is indeed continuous at 
L = 27r. 



5 Conclusion 

We have presented a new method allowing the computation of the Kramers rate prefac- 
tor in situations where the transition state undergoes a bifurcation. In contrast with the 
quadratic case, the prefactor is no longer independent of the noise intensity e to leading 
order, but diverges like Ce~", where a is equal to 1/4 times the number of vanishing 
eigenvalues. The constant C can in fact be computed in a full neighborhood of the bifur- 
cation point, and involves universal functions, depending only on the type of bifurcation. 
A similar non-Arrhenius behavior of the prefactor has been observed in irreversible sys- 
tems |MS96j . but there it has an entirely different origin, namely the development of a 
caustic singularity in the most probable exit path. 

Acknowledgments. We would like to thank Dan Stein for helpful advice, and for shar- 
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"Spectral Structures and Topological Methods in Mathematics" . 
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